Strongly correlated fermions with nonlinear energy dispersion and spontaneous 

generation of anisotropic phases 
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Using the bosonization approach we study fermionic systems with a nonlinear dispersion relation 
in dimension d > 2. We explicitly show how the band curvature gives rise to interaction terms in 
the bosonic version of the model. Although these terms are perturbatively irrelevant in relation 
to the Landau Fermi liquid fixed point, they become relevant perturbations when instabilities take 
place. Using a coherent state path integral technique we built up the effective action that gov- 
erns the dynamics of the Fermi surface fluctuations. We consider the combined effect of fermionic 
interactions and band curvature on possible anisotropic phases triggered by negative Landau pa- 
rameters (Pomeranchuck instabilities) . In particular we study in some detail the phase diagram for 
the isotropic/nematic/hexatic quantum phase transition. 
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I. INTRODUCTION 

One commonly used approximation to study strongly 
correlated fermions at low energies is the linearization 
of the fermion energy dispersion relation near the Fermi 
surface. For example, its implementation in the context 
of bosonization leads, together with the Renormalization 
Group (RG), to a powerful nonperturbative technique 
to deal with many-body problems. In one dimension 
this approximation gives rise to the Tomonaga-Luttinger 
model[l], while in d > 1 the Landau theory of Fermi liq- 
uids comes up as a fixed point in RG sense [2, 3, 4, 5]. In 
both cases, the bosonized Hamiltonian is quadratic and 
the model can be solved exactly, while small perturba- 
tions can be studied using RG. In particular, nonlinear 
terms in the dispersion relation are perturbatively irrel- 
evant, that is, they do not modify the long- wavelength 
properties of the system. 

However, in some cases, the linear approximation must 
be improved. Induced nonlinear terms in the quasipar- 
ticle energy dispersion become important when fermions 
are coupled to transverse fluctuating gauge fields. This 
is the case of some models of high Tc superconductors [6] 
and gauge theories of the half-filled {v — 1/2) Quantum 
Hall Effect [7]. 

More recently, it was suggested the possibility of hav- 
ing anisotropic ground states driven by spontaneous rota- 
tional symmetry breaking (Quantum Liquid Crystals) [8, 
9] . These new phases were proposed to describe transport 
properties of half- filled quantum Hall systems [9] and high 
Tc superconductors [8]. They can be associated to Pomer- 
anchuk instabilities [10] of the isotropic Fermi surface and 
for these novel ground states become stable a nonlinear 
fermion dispersion relation turns out to be essential[ll]. 
An interesting related phenomenon was also pointed out 



in ref. [12], where a Hubbard model is studied using 
RG techniques. There, it was shown that for a certain 
region of the parameter space strong forward scattering 
interactions favor Pomeranchuk instabilities leading to a 
breakdown of the discrete rotational symmetry. 

With these motivations in mind, we will present a sys- 
tematic study of the bosonization of fermionic systems 
with a nonlinear energy dispersion relation in any num- 
ber of dimensions. 

Concerning the Luttinger model case {d = 1), the first 
work were nonlinear dispersion terms were taken into ac- 
count was carried out by Haldane[l]. He showed that 
the resulting bosonized Hamiltonian is modified by the 
addition of nonquadratic terms in the bosonic variables. 
Explicit corrections to the corresponding one-particle 
Green function were recently computed in ref. [13]. For 
dimensions greater than one, the influence of nonlin- 
ear dispersion terms on the one-particle Green function 
was studied by Kopietz in the framework of functional 
bosonization[14] . 

In this work we are interested in studying the dynamics 
of the Fermi surface in dimension d > 2 and getting an 
explicit understanding of how the nonlinear dispersion 
relation could stabilize possible phases, other than the 
isotropic Fermi liquid one. 

The idea of a Fermi surface as a dynamical quantum 
extended object was originally introduced by Luther[15] 
and improved by Haldane[16, 17]. This concept was de- 
veloped in great detail by Castro Neto and Fradkin[2, 3] 
and by Houghton and Marston[4]. In ref. [2], the 
bosonized theory is written in a coherent state basis \(f> > 
representing deformations of the Fermi surface. In this 
way, the quantum dynamics of the system (the partition 
function) is expressed as a Feynman path integral where 
the "sum over paths" corresponds to summing up the 
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contributions coming from all possible deformations of 
the Fermi surface. 

The first part of this paper is devoted to apply the 
abovcmcntioned formalism to explicitly show how non- 
linear terms in the energy dispersion relation contribute 
with interacting nonquadratic terms in the bosonized ac- 
tion, in arbitrary dimensions. Eq. 4.7 below is one of 
the main results of this paper, showing the effective low 
energy Lagrangian of the system. There, (pq(k) repre- 
sents deformations of the Fermi surface at the Fermi 
point k (particle- hole excitations with momentum q). 
The second and third derivatives of the dispersion re- 
lation lead to the cubic and quartic bosonic terms, re- 
spectively. Notice that this general formulation stands 
for arbitrary smooth Fermi surfaces being particularly 
suitable for studying phases where shape deformations 
are present. 

These deformations can be classified according to the 
symmetries of an order parameter, similarly to the clas- 
sification of classical liquid erystals[18]. In fact, we can 
think about the quantum equivalent of smectic[19, 20], 
nematic[21] or hexatic[ll] phases and their corresponding 
quantum phase transitions. 

In ref. [11], V. Oganesian et al. showed for 
the first time that quantum isotropic/nematic and 
isotropic/hexatic phase transitions are possible in sys- 
tems where the Landau parameters F-2 and of the 
usual Fermi liquid theory[22] assume large and negative 
values (for a definition of the F„'s sec Eq. 5.5 below). 
In that reference, one important ingredient to stabilize 
the anisotropic states is the consideration of a nonlin- 
ear energy dispersion relation in the model Hamiltonian. 
The corresponding electronic properties are very promis- 
ing since the quantmn nematic and hcxatic states seem 
to present non- Fermi liquid behavior [11]. 

For these reasons, in the second part of this paper 
we use the nonperturbative bosonization approach to 
study quantum phase transitions to anisotropic electronic 
states in two dimensional systems. In particular, we will 
concentrate in nematic and hcxatic quantum liquid crys- 
tal phases where the order parameter is invariant under 
TT and 7r/3 rotations, respectively. 

The main result of this paper is displayed by the phase 
diagrams in figures 1 and 2. By integrating out all the 
stable modes we obtain an effective free energy at zero 
temperature as a function of the Landau parameters F2 
and Fq. We find different behaviors depending on the 
relative values of the stable Landau parameters F„ (n ^ 
2, 6). When these parameters are small, the phase dia- 
gram has a tricritical point where two second order phase 
transitions (isotropic/nematic, isotropic/hexatic) and a 
first order one (nematic/hexatic) meet together. How- 
ever if the stable F„'s are not small, a coexisting nematic- 
hexatic phase with a tetracritical point is possible. 

In the rest of the paper we explicitly develop the math- 
ematical details leading to these results. In section §11 
we present our Hamiltonian model for spinlcss fermions. 
In section §111 we show how to apply the bosonization 



method to Hamiltonians with a nonlinear dispersion re- 
lation and explicitly compute the corresponding non- 
quadratic bosonized terms. Then, in section §IV we ana- 
lyze the Fermi surface dynamics, building up a coherent- 
state path integral formulation for the partition function 
of the system. Finally, in section §V we analyze the possi- 
bility of isotropic-nematic-hexatic quantum phase transi- 
tions. Section §VI is devoted to a discussion of the results 
and to the presentation of our conclusions. 



II. THE HAMILTONIAN 

We consider a fermionic system characterized by a 
smooth Fermi surface given by the set of Fermi points 
kp satisfying e{kF) = H, where e{k) is an arbitrary en- 
ergy dispersion relation. The one particle excitations are 
associated with a set of operators ct and Cjj, creating 

and destroying a fermion with momentum k. These op- 
erators satisfy the usual fermionic anticommutation rela- 
tions. For simplicity we ignore the spin degree of freedom, 
however the extension to spinfull fermions is straightfor- 
ward. 

In general, the Hamiltonian can be written in the form 



H = Hq + i^int • 

The free (quadratic) term is given by 



(2.1) 



(2.2) 



A general two-body interaction term can be written as 



Hint Yj 4F,fe^('?i 4^ 



kF,k'„,q 



(2.3) 

where /^^ p (g) is the scattering amplitude among two 
particle-hole pairs with momentum q, at the Fermi points 
kp and k'p. 

In Eq. 2.2 the energy dispersion relation e{k) can be 
expanded in powers of q = k — kp, 
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3! dkidkjdki 



k=kF 



(2.4) 



where vp = vp{kp) = Ve(fc) 



^ ^ is the Fermi velocity. 

k=kF 



III. THE BOSONIZED HAMILTONIAN 

Bosonization is a powerful nonperturbative technique 
to deal with interacting fermions. In the case of two 



3 



dimensional parity breaking systems, it can be imple- 
mented in terms of a dual gauge theory (see for instance 
ref. [23] and references therein). On the other hand, 
the bosonization of a parity preserving system at finite 
density can be accomplished by introducing a restricted 
Hilbert space of small energy particle-hole fluctuations 
around the Fermi surface. In this case, the general for- 
malism was developed in refs. [2, 3, 4, 5]. 

In this section we find a bosonic representation for the 
Hamiltonian system 2.1 when a general energy dispersion 
relation (Eq. 2.4) is considered. In order to establish no- 
tation and to make this paper self-contained we will first 
summarize the main concepts of bosonization by follow- 
ing refs. [2, 3]. 

We define a reference state \FS) by applying fermionic 
creation operators to the vacuum state |0) so as to occupy 
all the states up to the Fermi surface, 



(3.1) 



We use this state to normal order all the relevant opera- 
tors of the theory according to 



: O -.= 0- {FS\d\FS). 



(3.2) 



The low energy behavior of the system is essentially 
described in terms of the particle-hole bosonic operator 



n,Kfc,t) = 4_,(f)cg^|(t) 



(3.3) 



where fc ~ fc^ and small q fluctuations are restricted to a 
thin shell around the Fermi surface. In fact, the approx- 
imation that defines the restricted Hilbert space of inter- 
est can be defined by the condition q < D < A « kp, 
where D is the shell thickness and A is the width of 
the finite amount of patches used to cover the Fermi 
surface[17]. These restrictions mean that the physical 
Hilbert space considered corresponds to a subset of ex- 
citations above \FS) mainly generated by small angle 
scattering processes. 

In this space the operators 3.3 satisfy the following 
commutation relation[2] 

[ng-(fc), {k')] = 5%^%,5q,q' q- vp 6 {i^ - eg) , (3.4) 

where because of the last delta function k is constrained 
to lie on the Fermi surface. For an arbitrary value of 
q, the operators n^{k,t) do not annihilate the reference 
state. However, we can define the operators. 



Hk,kF) 



t yjN{kF)V\q-VF\ 
\nS)0{q- n) + n-^{k)9{-q- n)} (3.5) 



4(^^) = E 



t yjN{kF)V\q-VF\ 
{n_^ik)eiq- n) + n+^{k)e{-q- n)} (3.6) 



(n is a unit vector normal to the Fermi surface at kp)- 
The smearing function $(A;, kp) is one, if k belongs to 
the patch labeled by kp and zero otherwise. In the ther- 
modynamic limit we have 



lim ^{k,kF) = Stt 



(3.7) 



the local density of states N{kF) is given by 

N{kF) = kF)fS{fi - e{k)), (3.8) 



and the operators a^{k), at(fc) satisfy 
a^kF)\FS}=0, 



(3.9) 



[aq{kF), ai, (kp)] = dj^^ p^ {S^^^ + 5^_^) , (3.10) 

generating the whole restricted Hilbert space of states. 
In this space, the fcrmion operator 

V(f,fcF) = ^e*«-%-(fcF) , q«k (3.11) 



can be written in bosonic form as[3] 



W,kF) = \ ^^U{kF)x 



a 



E 



-iq-r 



^ N{kF)V\q- VF 



■n-^{kF) 



,(3.12) 



or equivalently, in terms of a^{kF j, 



ip{f,kF) = 



'N(k, 



-Uikp) X 



E 

i3-n>0 



{e-'^^^ikp) - e^'^f'a^ikF)} 



\lN{kF)V\ 



q ■ Vf\ 



(3.13) 



where a is an ultraviolet cut-off and U{kF) are the 
"Klein Factors" that guarantee anticommutation rela- 
tions among operators with different kp (for an explicit 
expression of the Klein Factors see ref. [4]). 

Using the fermion-boson mapping 3.13 and the bosonic 
commutation relations 3.10 we can get the bosonized 
projection of any fermionic operator onto the restricted 
Hilbert space of states. For instance, the interacting part 
of the Hamiltonian (Eq. 2.3) is simply bosonized since it 
can be written in terms of nglkp), provided we restrict 
the momentum space considering only small angle scat- 
tering processes. Normal ordering the projected Eq. 2.3 



4 



we find, 

Hint = \Y. Pk.,k'^<l}\'i-M"'\<l-^F\"'^ 

+ a^(kF)a^{k'p) 9(-g • fcii-)6(g • k'p) + h.cj , 

(3.14) 

where we have introduced the adimensional Landau func- 
tion F^^^j:,^{q) = N'l\kp)Ny\k'p)t^^^^^.^{q), and 6 is 

the usual Heavisidc function. 

The restriction of the Hilbert space is usually justi- 
fied by RG arguments [24]. The long angle scattering 
coupling constants flow to zero at long distances, leav- 
ing a bosonized Hamiltonian only containing small angle 
scattering excitations. That is, the long angle scatter- 
ing operators arc perturbatively irrelevant in the Rcnor- 
malization Group sense. However, they rcnormalizc the 
parameters in the Hamiltonian. For this reason, the cou- 
plings Fj:^ p (q) in Eq. 3.14 should be considered as phe- 
nomenological inputs with no trivial connection with the 
microscopic ones. This limitation is at the heart of the 
bosonization procedure. However, this technique gives 
very general and powerful results concerning phase dia- 
grams and the universal structure of fermionic correlation 
functions. Of course, to make contact with microscopic 
models nontrivial numerical computations are necessary. 

Notice that for a fixed q, the first term in Eq. 3.14 
represents interactions among particle-hole pairs in the 
same hemisphere of the Fermi surface (with respect to 
the direction of q), while in the second term the interac- 
tion mixes the two hemispheres. This second term does 
not contribute to the asymptotic fermionic correlation 
functions [3], however, we will keep this term since it could 
become relevant in the case of nested Fermi surfaces. 

The bosonization of the free fermionic Hamiltonian is 
less trivial. Firstly, let us see what kind of terms appear 
in Hq (Eq. 2.2). The tensor structure of the Hamilto- 
nian is conveniently written in terms of a local reference 
frame defined by the unit vectors fi and t, normal and 
tangent to the Fermi surface, respectively. In terms of 
these directions it is easy to see that when considering 
the expansion 2.4, the first term in Hq only contains the 
normal field derivative 



(Tx ip'{x, kp) n • V ijj{x, kp). 



(3.15) 



For the second order derivatives of the dispersion relation 
we can write 



dkidkj 



^1 ninj+^2 {uitj +njti), (3.16) 



where ^i, ^2 and ^3 are functions of kp. A similar ex- 
pression can be written for the third order rank tensor 



containing the third order derivatives. Therefore, the 
X — space representation of the dispersive part of the 
Hamiltonian contains normal (fi • V) as well as tangential 
field derivatives (i • V). 

In order to obtain the bosonized form of Hq we consider 
the point-splitted product of the fermion field and its 
adjoint, along a general direction a. Using Eq. 3.12 and 
the Baker-Hausdorff formula 



(3.17) 



we obtain, 



1^ - 



kp^ (^f+ ia, kp^ := e^(«' x 



: e 



^ N{kp)Vq- VP 



sin ( '^Y~ni^{kp) 



where 



G{a, kp) = / 



^ (e*«-" - 1) 
N{^p)V\q-vp\ 



(3.18) 



(3.19) 



g-fi>0 



Let us first consider a direction a locally tangent to the 
Fermi surface (a = ef). In this case, when summing 
in Eq.3.19 over the normal {qj^ = q ■ h) and tangent 
{qr = q-t) components we get, 



G{e,kF) 



I — V(e-«--i) y 



N{kF)VvF % 



gw>0 



and considering eqx « 1, we obtain 



G{e,kp) = - 



2 rD 



(^A) . 

-1; / dqN -TZTT 

^vf Jo |9jv| 



(3.20) 



(3.21) 



(for an anisotropic Fermi surface Vp = \vp\ is kp de- 
pendent). Note that the remaining integral in q^ has a 
logarithmic infrared divergence, G(e, kp) oc — ln(y), and 
since the normal order in Eq.3.18 is a regular function 
we conclude that in the thermodynamic limit: 



22) 



This result implies that the tangent derivatives do not 
contribute to the projected bosonized Hamiltonian. In 
this regard, we recall that the bosonic particle-hole ex- 
citations tangent to the Fermi surface do not either 
contribute to the asymptotic form of the correlation 
function[3], however they do contribute to the density of 
states and they are crucial to obtain the correct specific 
heat and other thermodynamic properties [2]. 

We will now concentrate on the bosonization of 
fermionic terms containing normal derivatives. In this 
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case wc can generalize a calculation proposed by Hal- 
dane in one spatial dimension[l] to the case of arbitrary 
dimensions. Let us consider the integral, 

A = J df : — n, kp^ V (^+ n, ki 



(3.23) 

Introducing Eq. 3.11 into Eq. 3.23 and expanding in pow- 
ers of e we find on the one hand, 

(-0" 



On the other hand, choosing a = e n in Eq. 3.18, replac- 
ing into Eq. 3.23 and expanding in powers of e we find 
the bosonic version of A. Thus, comparing these two 
expressions order by order in e, we find the bosonized 
projected Hamiltonian. When expanding the dispersion 
relation up to third order in the derivatives we get 

HQ = Y,hi{kF) + h^{kF) + hSF), (3.25) 



A = Y^^-^e-Y.^q-hT ■.c\{kF)cSF):. (3.24) with 



hi 
h2 



VF 



\q-n\ ai{kF)ai^{kF) 



q-n>0 

13 



2{N{kF)VvFY' 



— J2 Wi ■ ^l'/^|92 • n|i/2|g3 • 5{qi +q2 + 9*3) x 



|at (fcF)at {kF)aq^ {kp) Q{-qi ■ n)Q{-q2 ■ n)e(g*3 • n) + h. c. | 



4! N{kF)VvF - 



|4 ai^{kF)a^2if^F)a^3{kF)a^i{kF)Q{-qi ■ n)Q{q2 ■ n)Q{q3 ■ h)Q{qi ■ h) 

+ 3 at (fcF)at ikF)a^s {kF)a^i (fc^) &{-qi ■ n)Q{-q2 ■ n)Q{q3 ■ n)&{q4, • n) + h. c. | , 



(3.26) 



(3.27) 



(3.28) 



where 



P = 



7 



d^ejk) 
dkjdki 

d^ejk) 
dkidkjdki 



k—kf 



niTijni 



(3.29) 
(3.30) 



k—kp 



The coefficient /? is related to particle-hole asymmetry. 
In fact, in a system invariant under charge conjugation, 
e{kF + q) — n is odd under the transformation q —q. 
In this case, the even derivatives of e at the Fermi level 
[k = kp) vanish. For stability reasons we will consider 
7 > 0. If 7 happens to be negative, then we should 
continue expanding the fermion dispersion relation until 
a well defined result be achieved. 

Eqs. 3.25 to 3.28 display the main results of this sec- 
tion. The fti-term corresponds to the bosonized free 
fermion Hamiltonian when a linear dispersion relation is 
considered and coincides with that one computed in rcfs. 
[2, 4]. Then, the final bosonized Hamiltonian contains 
a quadratic part (X^^t^- + Hint) plus a nonquadratic 
(/i2+/i3)-tcrm which is related to dispersive effects on free 
fermions. In order to calculate any observable the disper- 



sive part should be treated by means of perturbation the- 
ory. In the Fermi liquid regime the nonquadratic terms 
are irrelevant, although interesting effects were studied 
in the context of Landau theory [25]. 

Here, wc arc interested in understanding the role 
played by the presence of dispersion effects when the in- 
teracting system is otherwise unstable. The next two 
sections are devoted to study this issue. 



IV. DYNAMICS OF THE FERMI SURFACE 

As shown in ref. [2], the Fermi surface deformations 
can be associated to collective particle-hole excitations 
described by coherent sates in the bosonized theory. In 
the stable case, contact is made with the Landau theory 
for Fermi liquids. As we will sec, this procedure applies 
equally well to the general case where dispersion effects 
must be taken into account. 

Following ref. [2], wc define a many-body state which 
is a direct product of coherent states parametrized by a 
complex field (j}q{k): 



umFS), 



(4.1) 



6 



where 



and 



U{cl>) = e 



(4.2) 



1/2 



m ( — ] 

{(l)g-ikF)al{kF) - <l)}{kF)a^{kF)} . (4.3) 

The "deformation" fields satisfy (j)t(kF) = (l>-q{kF), im- 
plying r*(0) = -T{(l>) and C/t((/)) = t/-i ((/)). 

The system's partition function can be written in terms 
of this overcomplete coherent state basis by means of the 
path integral 



Z = j V(j)D(l)* 



(4.4) 



where S{(j),4'*) = J dt C{(p,(p* ) and 



t"^N{kF)V\q-VF\ 

kF,q 

{ct>\H\(^). 



mF,tf-^^ 



(4.5) 



The evaluation of {(j)\H\(j)) is straightforward, since the 
bosonized Hamiltonian (Eq. 3.25) is normal ordered and 
th{^ coherent states are eigenvalues of the destruction op- 
erator a^{kF)- 

UgikpMq- kF) + 'mpM-q- kF)) 
a,<kFM) = ^ , ^ 10) 



q-VF\ 



(4.6) 



Using Eqs. 3.14, 3.25 and 4.6, we find the fol- 
lowing Lagrangian (after the field redefinition (j) —>■ 

{VN{kF)vFY'^): 



qk 



dt 



+ 



qkk' 



+ ^ I 4 E "^fi (^)'^?3 + «2 + 93)+ 



3! 



kqi 



+ 5 E '^S-i (^^'^S-a (^)'^(^l + ... + 94) 



kqi 



(4.7) 



where k and k' lie on the Fermi surface and we have 
absorbed powers of N{Q)Vvf^&\ into the definition of /? 
and 7. 



The Lagrangian in Eq. 4.7 controls the low energy dy- 
namics of the Fermi surface. Note that this dynamics 
is extremely nonlocal, due to the factor l/|(f • n| in the 
kinetic term. In the region where the quadratic term 
is stable (positive definite), the first two lines represent 
the usual Landau theory of Fermi liquids. This is a 
fixed point in the RG sense, that is, the nonquadratic 
dispersive terms do not modify the asymptotic correla- 
tion functions. However, if interactions arc such that 
the quadratic term is unstable, the nonquadratic terms 
become relevant, stabilizing the theory in another fixed 
point. 

A word of caution is necessary to understand the valid- 
ity of Eq. 4.7. According to the Rcnormalization Group 
theory, not only the couplings Fj: j:, (q) but also the pa- 
rameters /3 and 7 will be renormalized by irrelevant op- 
erators. Therefore, although we showed how the band 
curvature generates the nonquadratic terms in Eq. 4.7 
(see Eqs. 3.29 and 3.30), the actual calculation of these 
parameters for a given band structure is not at all trivial. 

Also, in the present paper we are interested in under- 
standing how the band curvature can drive the system to 
a new rotational symmetry breaking gromid state. For 
this reason, in Eq. 4.7, we have only considered small 
angle scattering processes with small momentum trans- 
fer, disregarding any irrelevant term coming from the 
interaction Hamiltonian. On the other hand, it is well 
known that some interaction channels, although irrele- 
vant, could give raised to different instabilities. For in- 
stance, the Kohn-Luttinger instability[27] comes from the 
competition of forward scattering with a BCS channel, 
even for repulsive interactions. Recently, a similar dy- 
namical effect in 2d Fermi liquids was reported in ref. 
[28] . The interplay among the anisotropic phases studied 
here and other interaction channels is a very interesting 
issue, and the multidimensional bosonization technique 
described in this paper seems to be a promising tool to 
face it. 

With this comments in mind, we can study for instance 
the static deformations of the Fermi surface described by 
the system's free energy, which can be computed as the 
action per unit time, when setting to zero the kinetic term 
in Eq. 4.7. Introducing the Fourier transformed field 

<fi{x,k) = j -0yi Mk) e-'^-^ (4.8) 

at each point of the Fermi surface, the expression for the 
free energy is simplified to 

F = ^ j dxdx'Y^ip{x,k)M{k-k',x-x')ip{x',k') 

kk' 

+ E/ dx\^^^{x,kf + l^{x,kfY (4.9) 



where 



M{k-k',x-x') 
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FIG. 1: Phase diagram of an electron gas with nematic and 
hexatic instabilities. In this case \Fn \ « 1, for n ^2,Q. The 
bold and dot lines represent second and first order transitions, 
respectively. The origin BP is a bicritical point 

(4.10) 

In the next section, using the framework provided by 
bosonization, the problem of spontaneously broken rota- 
tional invariance is explored. 

V. THE ISOTROPIC-NEMATIC-HEXATIC 
QUANTUM PHASE TRANSITIONS 

The order parameter of the two dimensional nematic 

phase is a second-rank traceless antisymmetric tensor. It 
is odd under 7r/2 spatial rotations, and the two indepen- 
dent components can bo cast in the form of a "headless" 
two dimensional vector called "director" . The director is 
invariant under a tt rotation, characterizing in this way 
the nematic statc[18]. On the other hand, an hexatic 
phase is characterized by an order parameter that is in- 
variant under 7r/3 rotations. These states are of great 
interest as they arc candidates to present non-Fermi liq- 
uid behavior in two dimensions[ll]. 

By means of the bosonization approach, we will study 
them as particular instabilities of an isotropic Fermi sur- 
face. The important role played by the fermionic inter- 



actions as well as the dispersion effects for the establish- 
ment of these phases will become explicit. 

Let us begin by considering an isotropic Fermi surface 
in two dimensions and a short ranged interaction. In 
this case, the Landau function only depends on the angle 
between the two Fermi vectors k and k' 

Fj:j:,{x-x')=Fe_e' S{x-x'). (5.1) 

With these considerations, the ground state of Eq. 4.9 
becomes homogeneous, {ip{x,k)) = ip{6), and the free 
energy is simplified to 

Now, introducing the representation 

^{0) = ^<^„e»« (5.3) 

n 

F{0) = ^F„e-^ (5.4) 

n 

((/?* = <fi-n) and considering a particle-hole symmetric 
system (/3 = 0) we obtain, 

n nmlp 

(5.5) 

In the case where Fn > —1, it is clear that the only 
minimum of Eq. 5.5 corresponds to (pn = 0, for all n. This 
is the usual criterium for the Fermi liquid stability[2, 10, 
16, 29]. However, if for some channel n it happens that 
Fn < —1, then it is possible to develop quantum phase 
transitions to anisotropic states. Suppose for instance 
that we have a ground state characterized by (^2 7^ 0, 
and (fin = for all n 7^ 2. In this case, the deformation 
field would have the form 

^{0)=4<fi2 {2cos2(e)-l}, (5.6) 

characterizing a two dimensional nematic order [18]. This 
would correspond to a transition where, because of the 
interactions, the otherwise circular shape of the Fermi 
surface is deformed into an ellipse. 

In general, as we mentioned above, it is interesting to 
study instabilities in the nematic/hexatic sector. They 
are expected to occur when F2 and Fq are large and 
negative[ll]. This case can be treated by integrating 
over all the stable modes (pn, n ^ 2, 6. As a result 
of this integration all the relevant terms compatible with 
the symmetry will be generated, obtaining the following 
effective free energy for the modes ^2 and 

F = \2ipl + Kifl + ^ {^2 + '^t] + mflfl, (5-7) 
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FIG. 2: Phase diagram of an electron gas with nematic and 
hexatic instabilities. In this case, the Landau parameters are 

not restricted to small values and /i_R > The bold 

lines represent second order phase transitions and TB is a 
tetracritical point 



which is written in terms of renormalizcd interaction con- 
stants (A2 = + ^2), \ = vf{'^ + Fa). This free en- 
ergy presents the interesting possibility of quantum phase 
transitions [30] among isotropic (1^2 = ^f, = 0); nematic 
((^2 7^ , = 0) and hexatic states ((^2 = , 7^ 0). 
There is also the possibility of coexisting nematic-hexatic 
phases (</?2 7^ , (pQ^Q). 

It is important to note that while the phase transitions 
are triggered by the values of F2 and Fg, the structure of 
the phase diagram depends on the other Landau param- 
eters F„ (n 7^ 2, 6) hidden in and hr. For instance, 
if |F„| << 1 (for all n ^ 2,6), then /ijj « 7fl/2, and 
the qualitative phase diagram is that shown in figure 1. 
This diagram presents a bicritical point at the origin of 
the Ai — A2 plane. The bold lines represent second or- 
der phase transitions between the isotropic/nematic and 
the isotropic/hexatic phases. Moreover, the interphase 
nematic/hcxatic (doted line in Fig. 1) corresponds to a 
first order phase transition. 

However, if F„ « 1 there is the possibility of hav- 
ing 7^/16 < /i/j. In this case the phase diagram 
changes qualitatively (see Fig. 2). There is a tetracriti- 
cal point with four second order phase transitions. For 
Ai < 0, A2 < there is a region of coexisting phases 
which is absent in the preceding case. 



In order to study the electronic properties of these 
phases we need to consider the dynamics associated with 
the Lagrangian Eq. 4.7 and to evaluate quantum fluctu- 
ations around the saddle points found in these sections. 
This discussion will be presented elsewhere. 



VI. DISCUSSION AND CONCLUSIONS 

In this work we have constructed the bosonized ac- 
tion for a general fermionic system having a smooth 
Fermi surface and a nonlinear energy dispersion relation 
in higher dimensions. We have shown that the effect of 
the nonlinear terms in the energy dispersion is that of 
producing interactions in the bosonized theory. To the 
best of our knowledge this is the first explicit generaliza- 
tion to higher dimensions of the well known analogous 
result in one dimension[l]. 

The Hamiltonian 3.25, the corresponding action in the 
coherent state basis 4.7 and the free energy in Eq. 4.9 are 
the main results of this paper. 

From a physical point of view, dispersion effects are 
irrelevant when the system is in the normal Fermi liquid 
regime, however, they are essential when due to some 
fcrmion interactions an instability of the Fermi surface 
occms, driving the electronic system outside this regime. 
In this later case, the induced nonquadratic terms in the 
bosonized free energy 4.9 will stabilize the electronic sys- 
tem in a new ground state. 

In particular we have concentrated in two quantum liq- 
uid crystal states: nematic and hexatic. The correspond- 
ing phase transitions are triggered by negative values of 
the Landau parameters F2 and Fq, associated with the 
fermion interaction [11]. However, the qualitative struc- 
ture of the phase diagrams (Figs. 1 and 2) depends on 
the relative values of all the other "stable" Landau pa- 
rameters. When they are small, the phase diagram has 
a tricritical point where two second order phase tran- 
sitions (isotropic/nematic, isotropic/hexatic) and a first 
order one (nematic/hcxatic) meet together. In the oppo- 
site case, there is the possibility of having a coexisting 
nematic-hexatic phase with a tetracritical point. 

These results were obtained by means of a mean field 
saddle-point calculation on the bosonized action. It is 
clear that some of the global features (specially the first 
order phase transition in Fig. 1) could be modified by 
quantum fluctuations. 

In particular, it would be very interesting to investigate 
the dynamical electronic properties of these phases. In 
this regard, it was shown[ll] that the Goldstone modes of 
the nematic phase are damped (except for certain Fermi 
points dictated by symmetry considerations) while the 
fermion correlation function shows non- Fermi Liquid be- 
havior. In order to address these points in the nonper- 
turbative bosonization framework it is necessary to intro- 
duce quantum fluctuations around each saddle-point and 
employ the Lagrangian 4.7 and the fermionic operator 
3.13 to compute fermion correlators. We hope to present 
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